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Casimir interactions between macroscopic objects are strongly influenced by their 
geometrical features as shape and orientation as well as by their material properties. 
The effect of geometry is commonly obtained from the proximity approximation 
(PA). Here we present a path integral quantization for the electromagnetic field 
in the presence of deformed metallic surfaces. From the resulting effective action 
the Casimir force between the surfaces can be calculated without the PA. For 
corrugated surfaces the force is obtained both perturbatively for small deformations 
and by a numerical approach for general deformation amplitudes. For general 
dielectric materials with flat surfaces a path integral based derivation of the Lifshitz 
theory is outlined, pointing towards a possible approach to study the combined 
effect of deformations and material properties. 



1 Introduction and Summary 

Casimir's seminal prediction of an universal attractive force between two un- 
charged metallic plates has been confirmed recently in the distance range from 
lOOnm to several micromet er with increa sing accuracy in a number of beautiful 



high precision experiments UJaMcaasail High accuracy could be only achieved 
by a careful inclusion of the corrections to Casimir's ideal result from finite 
conductivity of the materials, surface roughness and non-zero temperature.^ 
Except one, ® the above cited experiments deviate from the flat-plate geome- 
try in that they use a plane-sphere configuration in order to avoid the difficulty 
of maintaining a constant plate separation over the entire surface area. The 
plane-sphere configuration is the standard situation where the proximity (or 
Derjaguin) approximation (PA) 1*3 is commonly employed to calculate geome- 
try induced changes of the force. Consider a non-planar surface above a planar 
surface where their local distance, measured normal to the planar surface, is 
given by H + /i(y||) with the 2D in-plane coordinates, i.e., h(y\\) — cor- 
responds to two parallel flat plates. The PA assumes that the Casimir energy 
can be computed as the sum of local contributions between flat surface ele- 
ments at their local distance, yielding the total energy per unit surface area 

£pa= / d 2 yil £ a [H + h( yil )}, (1) 

where £q(H) = — (ir 2 / 720) (He/ H 3 ) is the Casimir energy per unit area be- 
tween two plane surfaces. This approximation is ambiguous since one could 
just as well measure the local distances normal to the curved surface, yielding 
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a different result. Below, we will discuss this point further when compar- 
ing the PA to our perturbative and numerical results. The PA is expected 
to hold if both the local surface curvature radii are much larger than the 
local distance and strongly non-parallel surface elements are at larger separa- 
tions than more parallel onesM^J A conceptual different approximation is the 

1(51 

pairwisc summation of renormalizcd retarded van der Waals forces. For the 
geometry considered here it yields exactly the PA result of Eq. (Q.^ For the 
plane-sphere geometry used in experiments, one can set h{y\\) = yj/2R with 
radius of curvature R, yielding the approximate force Fpa = 2ttR£q(H). The 
latter formula is used in experiments, and there it is justified for that R is 
much larger than H. a 

For more general geometries the collective nature of the Casimir force 
suggests it to have a non-trivial and unexpected dependence on the shape of 
the interacting objects. Whereas the van der Waals force between electrically 
polarizable particles is always attractive, even the sign of the Casimir force 
is geometry dependent. There is little intuition for the value of the sign as 
demonstrated by the repulsive force for a thin conducting shell!™ Due to 
the importance of Casimir forces in many fundamental and applied contexts, 
it is highly desirable to obtain a better understanding of its strong geom- 
etry dependence, including distinctions from pairwise additive interactions 
and possibly repulsive forces between disconnected surfaces. A promising and 
experimentally testable route to this end is via modifications of the parallel 
plate geometry.El i n an experimental search for novel shape dependencies, 
Mohideen et al. measured the force between a sphere and a sinusoidally 
corrugated plate with a wav e length being larger than the studied range of 
sphere-plate separations.^! Their results showed clear deviations from the 
PA prediction. While it has been suggested that lateral forces caused this 
deviationJl^l, our results indicate a much stronger sensitivity to geometry at 
smaller corrugation wave lengths. 

We demonstrate that a path integral quantization of the electromagnetic 
field subject to approp riate bounda ry c onditions is a powerful method both 
for perturbativeli^^l and numerically computations of Casimir interactions. 
While the well-known boundary conditions for ideal metals are local, we derive 
non-local boundary conditions for general dielectric materials which reproduce 
the standard Lifshitz theory in the limit of flat surfaces. For the force between 
deformed surfaces of ideal metals we find the following results. For general 
uniaxial deformations we obtained the Casimir force perturbatively to second 
order in the deformation amplitude a. At this order, significant deviations 
from the PA result are found if the deformation wavelength A is of the same 
order or smaller than the mean surface distance H. We apply the perturba- 
tive results to sinusoidally corrugated surfaces, see Fig.^a). The deformation 



a A numerical scalar field analysis yields sizable deviations from the PA for R < 50 H 
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induced change in the force shows for A < If a slower decay ~ H~ 5 com- 
pared to the H~ 6 behavior from the PA. However, the perturbative approach 
is limited to the cases where a is smaller than both H and A. This can be 
understood by comparing the perturb ativ e approach to the multiple scatter- 
ing approach of Balian and Duplantiep^ 1 which yields the Casimir energy as 
a sum over waves propagating along closed paths which scatter at an even 
number of positions on the surfaces. Each path contributes with a geometric 
factor which measures at each position the scattering direction relative to the 
local surface normal. For a series of scatterings on the same surface with all 
scattering directions being almost perpendicular to the surface normal the 
geometric factor is small. Our perturbative approach allows in principle for 
paths with an arbitrarily large number of scatterings. However, it corresponds 
to an expansion of the geometric factor to second order in the surface defor- 
mations. This amounts to the restriction to two successive scatterings on one 
surface before the wave propagates to the other surface and experiences for 
the rest of its path only scatterings normal to the surface. This expansion 
is only justified for surface geometries where a is smaller than A so that the 
scattering directions lie nearly in the mean surface plane. The restriction 
can be bypassed by an exact evaluation of th e pa th integral by an algorithm 
which was developed for corrugated surfacesM^J This numerical approach is 
applied to the rectangular grating of Fig. ^b). It confirms the crossover of 
the force change from a decay ~ H~ 6 to ~ H~ 5 at H ~ A. In the limit of 
A < a strong corrections to the perturbative results are found. For A — > the 
force corresponds to that of two planar plates at a reduced distance H — a 
and provides an upper bound to the force at fixed H . For large A both the 
perturbative and the numerical results approach the PA prediction. However, 
we find that corrections to the PA decay slower with increasing A/a for the 
profile with edges [FigQIb)] as compared to the smooth profile [FigQJa)]. For 
this finding we provide an interpretation in terms of classical ray optics 




Figure 1. Geometry composed of a flat and a corrugated surface which is studied by per- 
turbation theory (a) and a numerical approach (b). For this two cases the corrections to 
the proximity approximation scale differently to zero for increasing A. 
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2 Path Integral Approach 



Due to the lack of space we give only a very brief account of the path integral 
approacrl^ and focus instead on explicit results for corrugated surfaces. For 
the uniaxial geometries considered here, the electromagnetic field can be split 
into transversal electric (TE) and magnetics (TM) modes which can both be 
described by a scalar field $. After a Wick rotation to imaginary time, the 
field is quantized using the Euclidean action 

S[$] = i /" d 4 X(V$) 2 . (2) 

The change in the ground state energy due to the presence of boundaries 
depends on the type of modes. TM modes are described by Dirichlet bound- 
ary conditions, <&\g — 0, and TE modes correspond to Neumann boundary 
conditions, dn'&ls — 0, where <9 n is the normal derivative on the surface point- 
ing into the space between the plates. The Casimir interaction between the 
two surfaces due to TM modes can be obtained from the restricted partition 
function 

Ztm = Z- 1 /»]][] 5MX a )}e- s ^\ (3) 

where Zq is the partition function in the absence of boundaries. For TE 
modes one obtains -Zte by replacing the argument of the delta-function 
by d n Q(X a ). The position of the surfaces are parameterized in the 4D 
Euclidean space by X x = [y,h(y\\)}, X 2 = [y,H], y = (ict = y ,y\\), 
where H is their mean distance, and the upper surface is planar. The 
partition functions can be computed by expanding the delta-functions as 
Y[ x ^[^K-^a)] = / T>ijj a exp[i J s dX a ip a $] with an auxiliary field tp a defined 
on each surface. Via the relation F = he/ (AL)dn InZ, with L the Euclidean 
system size in time direction, one finds for the force F per unit area the result 



F = -^ If Tr{M- 1 d H M). (4) 

This formula holds separately for both types of modes. For TM and TE 
modes, respectively, the 2x2 matrix kernel read 

M^ M (y,y') = G[X a {y)-X p {y')] (5a) 
^™(y,y') = d na{y{0 d npiA) G[X a {y) - X {y')\ (5b) 

with G(y,z) = (y 2 + z 2 ) -1 /47r 2 the free Green's function. Below, we will 
evaluate the force from Eq. |0J both perturbatively in /i(yy) and numerically. 

So far we applied the path integral approach only to ideal metal surfaces. 
If one assumes a splitting into TE and TM modes, arbitrary dielectric surfaces 
can be described by a path integral with non-local boundary conditions for the 
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scalar field^"". Planar surfaces with dielectric function t(w) are described by a 
mixed boundary condition which reads in momentum (fco,kii) representation 

[i-ra n ]$| s = o (6) 

with r = [e(i/c )fc2+k|]- 1 /2 and r = e (i/c )[ e (i/fc )fc2 +k 2]-i/2 for TM and TE 
modes, respectively. In the limit of ideal metals, e — + oo, this condition reduces 
to the above treated Dirichlct and Neumann boundary condition for TM and 
TE modes, respectively. The condition of Eq. 10 together with the action of 
Eq. (j2J provide a simple derivation^ of the standard Lifshitz theoryPH More 
interestingly, the application of more general non-local boundary conditions to 
gauge field path integrals might prove useful for studying correlations between 
geometrical and material dependencies of Casimir forces. 



3 Perturbative Results 



The logarithm of the partition function of Eq. J2J can be expanded in the 
height profileP^With the choice that / d 2 y\\ h(yn ) = there is no contribution 
at linear order in h(y»). At second order the iJ-dependent part of InZ has 
the form 

5{2)z = 2ib^ / d2y " h2{ y\\)-\ J An J d 2 y\ *(|yn-y|,|) [h(y l{ ) - My|,)] 2 

(7) 

with a kernel K which is different for TM and TE modesP^From this general 
result the force for the sinusoidal geometry of Fig. ^a) can be computed. 
Assuming h(y\\) = acos(2nyi/X), the total force from TM and TE modes 
assumes the form 



F = F Ra 



(8) 



with Ffl a t = — (n 2 / 240) he/ H 4 the force per unit area bet ween t wo flat plates. 
Although the full form of the function G(H/X) is available* I^ I IOJ we focus here 
on two limiting cases. For small and large corrugation lengths A, one finds 
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F/F aat = 1 + — — for A < H , F/F flat = 1 + 5-=j for A » H . (9) 



Therefore, for small A the change of the force decays only like H~ 
pared to ~ H~ 2 in the previously known limit of large A for which the PA is 
expected to hold. The divergence for small A is an artifact of the perturbative 
approach as we will show below. Next we compare these results to the pre- 
dictions of the PA. Depending on which surface this approximation is based, 
one obtains from Eq. to second order in a for the force 

2 2 2 

^PA.fW-Fflat = 1 + 5-^ , F P A,corr/F fla t = 1 + - 3tt 2 ^ (10) 
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for the flat and corrugated plate based PA, respectively. As expected, the 
later result is smaller due to the increased surface separations when measured 
normal to the corrugated plate. For A — > oo both results coincide and agree 
also with the perturbative finding. It is instructive to study how the difference 
between the full perturbative result of Eq. |JSJ and the PA decays to zero at 
large A. For this difference we find the following expressions which have 
opposite sign, 



F flat V 3 / A 2 ' F flat V 3 / A 

Thus, the flat surface based PA overestimates the force while the corrugated 
surface based PA yields a force which is too small. This underlines the ambigu- 
ity of the PA, even for large A or small surfac e cu rvature. Similar observations 
have been made for a sphere-plane geometryU^J However, even the exponent 
of a I A can depend on details of the corrugated surface. As we will show be- 
low, for the rectangular corrugation of Fig. ^b) the decay towards the PA 
result is slower as for the present geometry. In the next section we provide an 
explanation of this strong geometry sensitivity in terms of geometric optics. 



4 Exact Numerical Results 



We have seen that the perturbative approach is limited to smooth surface 

profiles where the deformation amplitude a sets the smallest of the geometrical 

length scales. The condition for the applicability of perturbation theory is 

that two arbitrary points on the deformed surface are connected by a vector 

which is nearly parallel to the reference plane. This is certainly not true 

for the profile of Fig. da) if A C a but is more generally violated if the 

profile has vertical segments as it is the case for the rectangular corrugation 

of Fig. ^b). For the latter profile, the perturbative expression of Eq. JJJ) 

indeed diverges. However, the general result of the path integral approach, 

Eq. can be used for a precise numerical computation of the force even 

in cases where perturbation theory fails. For periodic profiles one can apply 
[T"71 

the following approach. First, the Fourier transformed kernel Af(y,y') of 
Eq. JSJ is decomposed into contributions from momenta which are integer 
multiples of 2n/X, i.e., 

oo 

M(p,q)= N m (p ± , Pl )(2TT) 3 S(p ± +q ± )5( Pl +q 1 +2TTm/X) (12) 

m— — oo 

with = (po,p2). The 2x2 matrices N m are completely determined by a 
given uniaxial surface profile h(yi). Using Eq. Q), the force per unit area can 
be then obtained from 

he f°° r /x 

F = J dp ± p±_ J dpig(p±,pi), (13) 
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with 



g{p±,Pi) = tr (B 1 {p±,Pi)d H B(p ± ,p 1 )) , 



(14) 



where the lower case symbol tr denotes a partial trace with respect to the 
discrete indices k, I, and the matrix B is given by Bki(p±,pi) = Nk-i(p±,pi + 
2nl/X) with the N m defined by Eq. (|12() . In the following, we will apply this 
formula to the profile of Fig.^b) for which the matrices N m can be obtained 
analytically!"^ 

For A — ► 0, i.e., if A is the smallest length scale, the function g(p±,pi) can 
be calculated analytically, leading to g(p±,pi) = 2p/(e 2p ^ _a - ) — 1) for both 
TM and TE modes. After insertion into Eq. (|13|) . one finds the simple result 



He 



480 (H - a) 4 



(15) 



for the force from TM or TE modes, respectively. This result can be obviously 
interpreted as the force between two planar plates at a reduced distance H — a 
which reflects that the relevant modes do not " feel" the narrow valleys of the 
profile in the limit A <C a, H — a. In the opposite limit of very large A, the PA 
should be reliable, with the result 



F VK = - — - 
480 2 



1 



1 



(H + a) 4 {H 



(16) 



for both TM and TE modes. Here the PA yields the same result for the flat 
and the corrugated surface based approximation. 
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Figure 2. Casimir force -Ftm from TM modes for the geometry of Fig.^b). (a) Relative 
change compared to the force between two fiat plates for different corrugation lengths A. 
(b) Crossover between the proximity approximation at large A [Eq. 1161 1 and the reduced 
distance result of Eq. 1151 at small A. 
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For all intermediate values of A, the force can be calculated numerically 
from Eq. I|13|) by defining the M-th order approximation gAi(p±,Pi) to the 
function g(p_i_,Pi) by truncating the matrix Bki symmetrically around (k, I) — 
(0, 0) at order M so that the trace in Eq. ijTijl extends only over k, I = — (M — 
l)/2, . . . , (M — l)/2. This in turn, by numerical integration of Eq. I|13|). yields 
a series of forces which can be extrapolated to the actual force at M — > oo. 
The results for TM modes have been obtained in Ref. and are shown in 
Fig. The force -Ftm at fixed H/a is found to be monotonous in A/a and 
its minimal value is given by the PA result Fpa of Eq. Ijlfijl . At small A/a the 
exact result of Eq. I|15(l is recovered, and the divergence seen in perturbation 
theory, Eq. Jjyl, is removed. However, in agreement with perturbation theory 
for the sinusoidal profile, the relative change i*TM/-Fflat — 1 of the force shows 
a crossover at A ~ H from a H~ 2 decay ati/<Atoai? _1 decay at larger 
H. Fig. Hfb) shows how the limits of Eqs. (JTSJ, l|16|) ar e approached with 
decreasing or increasing A at fixed H/a. In the range of studied values for 
H/a the crossover between the two limits appears at A/a of order 10. As 
shown by the inset of Fig. |2{b) for H/a = 10, we find that the corrections to 
the results of Eqs. I|15fl. ifTBj l obey power laws, 

Fp — Ftm ^ A Ftm — Fpa ^ a 

for A/a <C 1 and A/a ^> 1, respectively. For TE modes we expect in the 
first case a decay ~ (A/a) 1 / 2 while in the latter case the exponent will be the 
same as for TM modes. Comparing to Eq. (|lll) we observe that the PA limit is 
approached with a different exponent for large A. This scaling behavior can be 
understood from classical ray optics. Following a recently proposed "optimal" 
PA for a scalar field-^, the Casimir energy is not obtained from the normal 
distance based on one of the surfaces [cf . Eq. 0J] but from the shortest surface- 
to-surface ray (of length £(y\\,z)) through given positions (y||, z) between the 
two surfaces of area A, 

Of course, this is still an approximation even for small surface curvature since 
the surfaces are treated as locally flat. However, we expect that this approx- 
imation yields the correct scaling behavior for large A/a. For the rectangular 
corrugation the lengths of the shortest paths can be computed analytically. 
For small a/H the dominant deviation from the standard PA comes from rays 
through positions located in equal sectors of almost triangular cross section 
for which the rays end exactly on the upper edges of the corrugated pro- 
file. For large A the volume of a sector becomes independent of A. In this 
case, a simple calculation yields, using Eq. I|18J) . the correction to the PA 



force, (F op t — FpA)/-Fflat = -g - ya/Ha/X, which reproduces our numerically 
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observed decay ~ a/A of Eq. l(T7|) . How can this be reconciled with our per- 
turbative result of Eq. for the sinusoidal profile? The crucial point is that 
the latter profile has a finite slope across the entire surface. Thus, the shortest 
length i deviates for all positions between the surfaces from the distance used 
in the standard PA. Using Eq. |JT5J), a model calculation for a corrugation 
with a piecewise constant slope indeed reproduces a decay ~ (a/A) 2 of the 
correction to the standard PA result. 

So far, we presented numerical results for TM modes only. For TE modes 
the same approach can be applied to the matrix of Eq. |J5J. The full electrody- 
namic Casimir force is then obtained as -Ftm + -Fte • Instead of presenting the 
results for TE modes independently, it is more useful to look at their relative 
contribution compared to TM modes. Fig. displays the ratio Ftm/Fte of 
the force from TM and TE modes for different corrugation lengths. Since the 
force is again bounded from below by Fpa and from above by Fq, one has 
Ftm/Fte — » 1 for both A/a — » oo and A/a — > 0. However, for intermediate 
corrugation lengths the ratio can differ from one, with the amount depending 
on the mean surface distance H. We observe the tendency that TM mode 
contributions dominate at smaller A while for larger A the TE modes provide 
a larger force. In contrast, for a sinusoidal profile the force from TM modes 
is larger at all A, at least to second order in a. For very small separations 
H > a the segments of the rectangular corrugation which are closer to the flat 
surface yield the main contribution and the PA is expected to hold, implying 
Ftm/Fte ~ 1 which indeed can be observed in Fig. [31 In the opposite case of 
large H ^> a, one obtains the flat plate geometry, and again Ftm/Fte — > 1. 
The latter limit is approached for both types of corrugations with a power law 
~ (if/a) -1 which we find in perturbation theory as well as in the numerical 
approach. 
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